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New experiments were proposed recently to investigate the regime of cold atomic and molecular 
ion-atom collision processes in a special hybrid neutral-atom-ion trap under high vacuum conditions. 
The collisional cooling of laser pre-cooled Ca+ ions by ultracold Na atoms is being studied. Modeling 
this process requires knowledge of the radiative lifetime of the excited singlet A^E"*" state of the 
(NaCa)+ molecular system. We calculate the rate coefficient for radiative charge transfer using 
a semiclassical approach. The dipole radial matrix elements between the ground and the excited 
states, and the potential curves were calculated using Complete Active Space Self-Consistent field 
and MoUer-Plesset second order perturbation theory (CASSCF/MP2) with an extended Gaussian 
basis, 6-311-|-G(3df). The semiclassical charge transfer rate coefficient was averaged over a thermal 
Maxwellian distribution. In addition we also present elastic collision cross sections and the spin- 
exchange cross section. The rate coefficient for charge transfer was found to be 2.3 x 10~^^ cm^/sec, 
while those for the elastic and spin-exchange cross sections were found to be several orders of 
magnitude higher (1.1 x 10~* cm^/sec and 2.3 x 10~^ cm'^/sec, respectively). This confirms our 
assumption that the milli-Kelvin regime of collisional cooling of calcium ions by sodium atoms is 
favorable with the respect to low loss of calcium ions due to the charge transfer. 

PACS numbers: 32.80.Pj, 34.80.Qb, 82.20.Pm, 82.30.Fi 



I. INTRODUCTION 



A broad range of techniques from atomic physics and 
optics have allowed the accurate manipulation of ul- 
tracold samples. Recent studies are probing ultracold 
atomic systems in which electric charges may play an 
important role 0, Q, Q. These include ultracold plas- 
mas 0, IE IE ) ultracold Rydber g ga ses H, , as well as 
ionization experiments in a BEC |1C| . 

Recently, we proposed experiments on simultaneously 
cooling and trapping ions and atoms in the same con- 
fined space The proposed experiments will investi- 
gate the regime of cold atomic and molecular ion-atom 
collision processes in a special hybrid neutral-atom-ion 
trap under high vacuum collisions. The idea is to use 
the techniques of laser cooling, trapping and manipula- 
tion of cold neutral atoms to make a refrigerator for the 
sympathetic cooling of an overlapping sample of atomic 
or molecular ions to very low temperatures. We intend 
to probe fundamental collision processes and cooling ki- 
netics in as low a temperature range as possible. 

The hybrid trap we have built consists of a conven- 
tional magneto-optic neutral atom trap (MOT) combined 
with a linear Paul radiofrequency quadruple trap. The 
linear Paul trap, with longitudinal electrostatic confine- 
ment, is chosen so as to suppress most of the micromotion 
and r.f. heating associated with conventional Paul traps 
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II. POTENTIAL CURVES AND DIPOLE 
MOMENT 

To calculate the cross sections and rate coefficients for 
the charge transfer and the elastic and spin-exchange col- 
lisions we have performed ab initio calculations of the 
adiabatic potential curves for (NaCa)+ molecular sys- 
tem in the X^S+, a^E+ and A^E+ states, which go 
asymptotically to the ground (Na+-I-Ca) and excited 
(Ca'*"-|-Na(3s)) states of the ion-atom quasimolecular sys- 
tem. The method used was second order Moller-Plesset 
perturbation theory (MP2) using a Gaussian triple 
zeta-|-diffuse-|-polarization basis set, 6-311+G(3df). This 
basis set included three sets of d-polarization functions 
and one set of /-polarization fimctions in addition to 
a diffuse function on each atomic center. The first ex- 
cited A^I]+ state was optimized for the second root of a 
CASSCF (Complete Active Space Self- Consistent Field) 
using eight orbitals in the active space. This was followed 
with an MP2 (Moller-Plesset second order perturbation 
theory) correction to the optimized CASSCF energy. 

Figure n shows the potential curves of the ground and 
excited states of (NaCa)+ molecular system in atomic 
units (a.u.) (see also Table P). The curves were obtained 
by performing a cubic spline fit to the calculated points. 
At small internuclear separations R, the potential curves 
were approximated with the short-range potential of the 
form V{R) = {A/ R) exp{-BR), where A and B are fit- 
ting parameters, and for R > 30 oq, the potential curves 
were fitted with the appropriate long-range potentials. 
For the excited states, the long-range potential was of 
the form 
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FIG. 1: Adiabatic potential curves for selected states of the 
(NaCa)"*" molecular system. Rg is a classical turning point. 
Shown is a free-bound radiative transition with the emission 
of a photon hv. 
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TABLE I: Potential energy and dipole matrix element val- 
ues obtained by cubic spline fit of the adiabatic potential 
and dipole matrix curves for selected states of the (NaCa)"*" 
molecular system. Potential curves are located relative to the 
ground asymptote. Notation 1.23[-4] stands for 1.23 x 10"'' 
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3.54[-l] 


4.16[-1] 


1.25 
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7.42 [-2] 


1.98[-1] 
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6.88[-2] 




5 


-2.79[-2] 
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dipole radial matrix elements at smaller distances. At 
small distances there is little overlap between the wave- 
functions of the ground and excited states for the ener- 
gies of the approaching sodium atoms considered here. 
At milli-Kelvin energies of the sodium atoms, the classi- 
cal turning point Rq was found to be always to the right 
of i? = 8 ao. 



III. CROSS SECTIONS 



FIG. 2: Calculated radial transition dipole matrix elements. 
The matrix elements are calculated at distances R > S ao 
only, because the current problem did not require the knowl- 
edge of the matrix elements at lower inter-nuclear distances 
due to small energies of the approaching sodium atoms. 

where the Na dipole (162.7 e^ag), quadrupole (1902 
e^flg), and octupole (55,518 e^al) polarizabihties were 
taken from Ref.0]. For the ground state X-'^E+, the form 
of the long-range potential was 

nm--l^. (2) 

where a = 156 e^ajj is the dipole polarizability of the 
calcium atom. 

Figure 121 shows the calculated radial dipole matrix el- 
ement between X-'^E"'" and A^E"*" states fitted with cubic 
splines (see also Table IJ) . The dipole matrix elements 
were calculated for R > 8 ao only, since for our purposes 
the problem did not require the exact knowledge of the 



A neutral atom of Na collides elastically with Ca+ in 
the singlet A^E^ or triplet a,^T,+ states of NaCa+. The 
corresponding elastic cross sections are given by 

(=0 

where k = \/2iiE/h with fi the reduced mass and E the 

S T 

collision energy in the center of mass system, and rj^ ' 
is the l^^ partial wave phase shift corresponding to the 
A^E+ and a'^E+ states, respectively. Another possible 
outcome of the collision between Na and Ca"*" is spin- 
exchange. In the elastic approximation lll^ . also known 
as the Degenerate Internal States (DIS) approximation 
[T^ , one describes this scattering process in terms of the 
singlet and triplet scattering phase shifts, and the spin- 
exchange cross section is given by 

oo 
i=0 
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Although this approximation is valid for collision energies 
larger than the hyperfine splitting, it usually gives the 
right order of magnitude at smaller energies |l5l | . 

S T" 

The phase shifts rj^ ' are determined from the contin- 
uum eigenfunctions ^^'^(i?), which are the regular solu- 
tions of the partial wave equation 



2fi 
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0. (5) 



The asymptotic form of y^''^{R) at large distances gives 
the elastic scattering phase shifts 



sm 



(6) 



The results of our calculations are shown on Figs. El 
and 2] for the singlet and triplet elastic cross sections, 
respectively, and on Fig. [3 for the spin-exchange cross 
section. In all cases, the s-wave contribution is dominant 
at energies corresponding to temperatures below 1 nK, 
in sharp contrast to neutral alkali atoms where this take 
place for energies around 100 /xK This is due to 

the very long range of the polarization potential as com- 
pared to the shorter range of van der Waals interactions 
between neutral atoms. Our ab initio potential curves 
are not of sufficient accuracy enough to predict the s- 
wave scattering lengths associated with the singlet and 
triplet potentials with a fair degree of certainty: the mea- 
surements that will be realized in our proposed hybrid 
MOT would help fine-tuning the potentials so that accu- 
rate estimates of the scattering lengths could be given. 
As E increases, many more partial waves contribute to 
the cross sections: already at E/ks ^ 1 /iK, 7 partial 
waves are necessary, and over 30 are required at 1 mK. 
Although shape resonances may play an important role 
at very small energies, their effect is not dominant for 
E/ks larger than 100 /j,K (see the smooth behavior of 
the cross sections in Figs|31-|n|. This rapid growth in 
the number of partial waves led us to consider semiclas- 
sical approximations for the cross sections. Following the 
method described in Q , we found that the cross sections 
are well represented by power-law functions of energy for 
the temperature range down to /iK. Figures 13 - [S] show 
the calculated collision cross sections together with the 
power-law fits 



<7cxch.{E) = Ccxch.E^^^"^ , 



(7) 
(8) 



where Cd. = 5310 and 5070 for the singlet A^J:+ and 
triplet a'^S^ states, respectively, and Ccxch. — 44. These 
coefficients were obtained by performing a linear fit of 
the cross sections plotted on a log-log scale in the energy 
range above 1 /xK. All quantities are expressed in atomic 
units. 

Using these approximate cross sections, we can cal- 
culate the various rate coefhcients TZ = (crv), where 
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FIG. 3: Elastic collision cross section for the A^E^ state of 
molecular system (NaCa)"''. The elastic cross section above 
10 fiK fits the approximate expression (j{E) = 5310-B~^''"^, 
where the energy is in units of Hartrees and the cross section 
is in units of Oq. 



V — \J2E I yu is the relative velocity, and (•••) implies 
averaging over the velocity distribution. Assuming a 
Maxwellian velocity distribution characterized by the 
temperature T, we obtain the following rate coefficients 



^cxcli. — 



r(^)(fcBT)i/6Cei., (9) 

TT/X V3/ 



cxch. 1 



(10) 



in atomic units. Note that while 7?.cxch. is independent 
of temperature, 'R-^x. varies very slowly with T. Figs. 
l^lshow that the power law fits are valid at temperatures 
1-10 /xK and higher, which validates the semiclassical 
approach to calculating the rate coefficients for elastic 
and spin-exchange collisions. 

As an example, suitable for the proposed experiment, 
we evaluated the rates for a temperature of 1 mK. They 
are T^-d. = 1.7 — 1.8 a.u. for elastic collisions and 7?,oxch. = 
0.38 a.u. for spin-exchange collisions. Converting to 
units cm'^ sec~^, we obtain TZ^x. = (1.05 — 1.10) x 10~^ 
cm'^ sec~^ and T^oxch. — 2.33 x 10^® cm"^ sec^^. For 
a typical MOT particle density n — 10^° cm"-^ these 
rates give the following values for effective elastic life- 
times 1/r: 1/tci. = 7^ol.n = (1.05 - 1.10) x 10^ sec"!, 
and 1/rcxch. = T^oxch.'^ = 23.3 sec~^ 



IV. RADIATIVE CHARGE TRANSFER 

The rate coefficient for radiative charge transfer may 
be expressed as an integral of the transition probability 
over the collision path, avera ged over the initial velocities 
and the impact parameters p?7| . Because many partial 
waves contribute even at 1 ^K, we use a semiclassical 
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FIG. 4: Elastic collision cross section for the a'^E"'" state of 
molecular system (NaCa)"*". The elastic cross section above 
10 fiK fits the approximate expression a{E) = 5070i5~^''"^, 
where the energy is in units of Hartrees and the cross section 
is in units of Oq. 
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FIG. 5: Spin-exchange collision cross section between states 
A^E"*" and a^S+ of the molecular system (NaCa) + . The cross 
section fits the expression <j{E) — 44i?~^/^, where the energy 
is in units of Hartrees and the cross section is in units of ag . 



treatment where the cross section for radiative charge 
transfer may be expressed as 



a{E) = 2tt 



1/2 



AiR) 



where b is the impact parameter, Rq is the classical turn- 
ing point (distance of closest approach), A{R) is the Ein- 
stein spontaneous emission transition probability, and 
V{R) is the entrance channel potential curve. Again, 
we can obtain a rate coefficient TZ^^ = (av) by averaging 
over a Maxwellian velocity distribution characterized by 



7^tr 80F 



1 
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dE db 
ksTJ Jo Jo 

bE^/^A{R)e-^/''^^ 



Iro '^%-byR^-ViR)/E]^/^- 
The transition probabilities are approximated by 

A{R) = ra^uj^{R)D{R) a.u.. 



(12) 



(13) 



where a is the fine structure constant, D{R) is the dipole 
moment (see Fig. \^ and 



hiu{R)^Vn{R)-Vn'{R) 



(14) 



is the energy difference between the entrance and exit 
potential curves (here A^E+ and X^E+, respectively), 
all at a separation R. Note that the effect of possible 
shape resonances is not taken into account in Ea. (|12|l : 
as mentioned in the previous section, shape resonances 
are not relevant for E/ks ~ fOO /iK or higher (although 
they may play a role at much lower energies). The in- 
tegrations of the Eq. I|12|l were evaluated numerically 
for temperatures ranging from 1 mK to 1000 K. The 
rate coefficient was found to be independent of tem- 
perature and was practically constant with the value of 
TZtv = 2.3 X 10~^^ cm'^/sec. With the attainable densities 
of ultracold sodium atoms of the order of n = 10^ — 10^^ 
cm'^, the rate coefficient gives for the the radiative life- 
time of the excited A^S+ state 



l/(7^t 



4 X 10" - 4 X 10^ 



(15) 



or in the range from several hours to several days. Paren- 
thetically, we note that, although charge transfer without 
photon emission is possible at higher collision energies 
(many eV), in the IK range and below, the small cou- 
pling between the a and X states of (NaCa)"*" (see Figure 
^ makes the rate for nonradiative charge transfer negli- 
gible. 

We also considered the possibility of a charge trans- 
fer process that could occur by laser excitation from the 
Pi}Ti^ state of NaCa"*" to a higher molecular state that 
correlates with Na+ (3s) plus an excited Ca level. The two 
wavelengths that must be considered are Na(3s — *■ 3p) at 
589.0 nm and Ca+(4s 4p) at 397.0 nm. An exami- 
nation of the asymptotic limits of molecular states lying 
above Na(3s)-|- Ca"'"(4s) indicates that the closest excited 
level would correspond to Na''"(3s) -I- Ca[^P(4p)] at 633.0 
nm, too far from resonance to be of importance. |20j| 



V. CONCLUSION 

In conclusion, we have shown that the semiclassical 
approach to evaluating the rate coefficients is a good 
approximation because many partial waves are involved 
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even at ^ 1 mK. The exeited singlet state of the col- 
liding molecular system (NaCa)+ is metastable against 
radiative charge transfer with the lifetime of the order 
of many hours. At the same time, the rates for elastic 
collisions and also for spin exchange are much greater, 
which provides an efficient mechanism for sympathetic 
collisional cooling of calcium ions by sodium atoms in 
the meV regime and below. We are setting up an exper- 
iment to measure collisional cooling by this mechanism. 
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